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Abstract 

We investigate the thermodynamics of Lovelock-Lifshitz black branes. We begin by introducing 
the finite action of third order Lovelock gravity in the presence of a massive vector field for a 
flat boundary, and use it to compute the energy density of these black branes. Using the field 
equations, we find a conserved quantity along the r coordinate that relates the metric parameters 
at the horizon and at infinity. Remarkably, though the subleading large-r behavior of Lovelock- 
Lifshitz black branes differs substantively from their Einsteinian Lifshitz counterparts, we find that 
the relationship between the energy density, temperature, and entropy density is unchanged from 
Einsteinian gravity. Using the first law of thermodynamics to obtain the relationship between 
entropy and temperature, we find that it too is the same as the Einsteinian case, apart from a 
constant of integration that depends on the Lovelock coefficients. 
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I. INTRODUCTION 



Gauge/gravity duality presents a powerful new framework with which one may study 
strongly coupled gauge theories. Indeed, it has been conjectured that there is a duality 
between certain strongly coupled gauge theories and weakly coupled string theories, which 
means that both theories describe the same physics, but calculations may be easier in one 
theory than the other. More specifically, the AdS/CFT correspondence relates A/" = 4 
supersymmetric SU(N) gauge theory to superstrings in 10 dimensions One typical 

and fascinating aspect of gauge/gravity duality is the property of holography, which states 
that the amount of information contained in the boundary gauge theory is the same as 
the one contained in the bulk string theory. This holographic toolbox has since spread its 
usage to many branches of physics from hydrodynamic to condensed matter systems 
^. However these studies have been concerned with relativistic theories with a conformal 
symmetry in the ultraviolet, which are described by asymptotically Anti-de Sitter (AdS) 
spacetimes. Applying gauge/gravity duality to condensed matter theories with anisotropic 
scaling symmetry, 

t XH, X ^ Ax (1) 

where z is the dynamical exponent, suggests that we can further extend holography to 
encompass nonrelativistic and Lifshitz field theories. 

From a holographic point of view, this suggests the following (asymptotic) form for the 
spacetime metric 



ds^ = (^-r^'df + ^ + r^dx^^ (2) 
that obeys the scale invariance 

t XH, r A"V, X Ax. (3) 

noted previously in a braneworld context ^ . A four-dimensional anisotropic scale invariant 
background using an action involving a two form and a three form field with a Chern-Simons 
coupling or a massive vector field can be engineered to yield solutions with this asymptotic 



behavior jc]. For these matter fields, lots of effort has been expended in extending this 
solution to the case of asymptotic Lifshitz solutions. One of the first analytic examples 
was reported in Ref. [tJ for a sort of higher- dimensional dilaton gravity without restricting 
the value of the dynamical exponent z. An exact topological black hole solution that is 



asymptotically Lifsliitz with z = 2 was obtained in jsj; further solutions with z = 4 and 
with spherical topology were subsequently obtained 9|. In general, however, such asymptotic 
Lifshitz black holes must be obtained numerically 8NlO|. Asymptotic Lifshitz solutions in 
the vacuum of higher-derivative gravity theories (with curvature-squared terms in the action) 
have been investigated ll| ; the higher-curvature terms with suitable coupling constant play 
the role of the desired matter. Recently, we introduced some solutions - both analytically 
and numerically - which can be regarded as higher-curvature modifications from Lovelock 
gravity to those obtained from Einsteinian gravity coupled to matter 12 1. 

The holography of gravity theories including higher powers of the curvature, particularly 
curvature-cubed interactions, have attracted increased attention 3, ll^- This is because, 
in the context of AdS / CFT correspondence, corrections from higher powers of the curvature 
must be considered on the gravity side of the correspondence in order to investigate CFTs 
with different values of their central charges. On the other side, the asymptotic Lifshitz 
black hole from a holographic point of view needs an extension of holography to encompass 
nonrelativistic field theories. Our aim here is to further develop the holographic dictionary 
for asymptotically Lifshitz spacetimes in Lovelock gravity by introducing the counterterm 
method for Lovelock-Lifshitz black branes 12| with fiat boundary. Indeed, the generalization 
of holographic techniques to this new context may offer new insights into the nature of the 
relation between quantum gravity in asymptotically non-AdS spacetimes and the dual field 
theory. Here we generalize the counterterm method introduced in 15| for four-dimensional 
Einstein gravity to (n-l- 1) -dimensional Lovelock gravity, which of course contains the higher- 
dimensional Einstein gravity and Gauss-Bonnet gravity too. The counterterm method for 
asymptotic AdS solutions of Lovelock gravity with flat boundary was introduced by us in 



16|. We generalize it to include asymptotic Lifshitz solutions of Lovelock gravity. 



We make use of these generalizations to investigate the thermodynamics of Lovelock- 
Lifshitz black branes [l2|. Somewhat remarkably we find that the relationship between the 
energy density, temperature, and entropy density is unchanged from Einsteinian gravity. We 
make use of the first law of thermodynamics to obtain the relationship between entropy and 
temperature. We find this is also the same as the Einsteinian case, apart from a constant 
of integration that depends on the parameters in the Lovelock action. 

The outline of our paper is as follows. We introduce the one dimensional Lagrangian of 
the Lovelock-Lifshitz black branes in Sec. HIl We also write down the field equations and 
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present the constant Cq which is preserved along the coordinate r. This constant relates the 
information at the horizon and at infinity. In Sec. IIIIl we generalize the counterterm method 



introduced for the four-dimensional Einstein equation [15| to Lovelock gravity, and present 
the finite action for Lovelock-Lifshitz black branes. Section |IV] is devoted to the calculation 
of the constant Cq at the horizon and at infinity. In Sec. El we obtain the energy density of 
the Lovelock-Lifshitz black branes in terms of the temperature and entropy and find that it 
is the same in Einstein [itI] and Lovelock gravities. We also consider the thermodynamics 
of the solutions. We finish our paper with some concluding remarks. 

II. FIELD EQUATIONS 

I — I 

The bulk action of third order Lovelock gravity [18] in the presence of an Abelian massive 
vector field A^^ may be written as 



/bulk = f rf"+^Xv^ (Cg + Lr, 

167r ./x 



Lg = £i + a;2£2 + as-^s - 2A, 

i:^m = -\f^uF^^ -]^m'A^A^^ (4) 

where F^y = d[^A^], A is the cosmological constant, 02 and 03 are second and third or- 
der Lovelock coefficients, £1 = i? is the Einstein-Hilbert Lagrangian, £2 = R^iv^sR^"'^'^ — 
AR^yR^^ + R^ is the second order Lovelock (Gauss- Bonnet) Lagrangian, and 



d3 I o TjiiuaK p TDpT I Q TDliv DCTK JDpT I O/i ryUfaK p pp 

3 — -ft + Zn"^ KaKpr^ + apJ^ ut^ fiK + '^^-ti J^aKupJ^ 

^3^^M.aK^^^^^ + 24:R^''"'''R^^R^, + im^^R.^R^ - URR^^'R^, (5) 

is the third order Lovelock Lagrangian. We assume that the Gauss-Bonnet coefficient, which 
has the dimension of (length)^, is positive as in heterotic string theory [l^. In Lovelock 
gravity only terms with order less than [{n + l)/2] (where [x] is the integer part of x) 
contribute to the field equations, the rest being total derivatives in the action. For third 
order Lovelock gravity we therefore consider {n + l)-dimensional spacetimes with n > 6 
(though in situations where we set ds = our solutions will be valid for n > 4). 

We write the spherically symmetric gauge field and metric of an (n -|- l)-dimensional 
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asymptotically Lifshitz static spacetime with zero curvature boundary as 

A = qe^^'^dt = q^—h{r)dt, 

n-l 



i=l 



S\2 



(6) 
(7) 

(8) 



i=l 



where the new metric functions are related to our previous notation 121] through the following 
relations 

1 T 

F{r) = -ln/(r) + 

1 r 
G{r) = — \n g{r) — In 



R{r] 



2 



H{r) = In h{r) + z In 



I 



(9) 



and in Eq. ([8]) we have chosen coordinates by taking R(r 



ln(r//). 



Our previous study of Lovelock-Lifshitz black holes 12| showed that if 
(n-l)z ^ 2(z-l)L^ 



m 



P 



A 



[{z - 1)2 + n{z - 2) + n^]L^ + n{n - l){a2f - 263; 



2/6 



(10) 



where = — 20:2/2 + 303, then the action (j4]) supports solutions asymptotic to the Lifshitz 
solution [l2| 

(11) 



1=1 



Since we are assuming spherical symmetry, we can reduce the action to one dimension and 
subsequently obtain the equations of motion. After integration by parts, the one dimensional 
Lagrangian may be written as Cid = + >Cim), where 

A 



£1, = {n-l){- 2^e2G + [2F'E' + {n - 2)E!^] - ^ [AF'B!^ + {n - A)B!^] e'^^^ 

J 
Clm = ig2(m2 + if'2g-2G)g-F+G+(n-l)K+2H^ 



(12) 
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where prime denotes the derivative with respect to r, and we define a2 = {n — 2)(n — 3)a2 
and ds = {n — 2)...{n — 5)0:3 for convenience. The equations of motion following from this 
action may be written as: 



Ir 



3 5 



F'+{n- 2)R' - -a2 [SF'R'^ + {n - A)R'^) e'^^ 
3 



/5\ ^-4G 



da {5F'R"^ + in- 6)R'^) e 



^F-G+{n-l)R 



= [2F'R' + {n- 2)R'^] - [AF' R!^ + {n - 4)i?'^] e 



(13) 

(14) 

(15) 
(16) 



2G 



+d3 [QF'R"' + {n- Q)R!^] e-^^ + 



2A 
n — 1 



,2G 



2(n - 



1 ^ ^ 



(17) 



Now subtracting the summation of Eqs. ( |T3l) and ( |T6i) from Eq. ( |T5l) one obtains: 

{2{F' - R') (1 - 2d2i?"e-2^ + 3«3i?"e-^^) e^-^+("-i)^ - g2^'e-^-^+("-i)^+2H|' ^ 
which shows that 

Co = 2(F' - R') (1 - 2d2i?''e-2« + Sa^R'^e-^"") e^-^+("-i)« - ^2^/g-F-G+(n-i)i?+2/f 



1 - 2^^? + 3^/ ) [rf + 2iz - l)f] - q\zh + Th')h 



,«3 2 



/4 



,n+2— 1 



is conserved along the radial coordinate r. One may note that the action ( 1121) and constant 
( IT8|) reduce to those of Ref. js] for n = 3 and 02 = 03 = 0. This conserved quantity is 
associated with the shift 

/F(r)\ 



R{r\ 
G{r. 
\H{r)J 



— ^ 



/ F{r) + 5\ 

G{r) 
\ H{r) + 5 J 



(19) 



where 5 is a constant, and reduces to the diffeomorphism introduced in |9!] for = 3 and 
"2 = 63 = 0. 

We pause to remark that for 2; = 1 with /(r) = g{r-)^ the constant ( IT8|) reduces to 



which is known to be constant in third order Lovelock gravity and is proportional to the 
mass parameter of the spacetime 20|. 



III. FINITE ACTION FOR LOVELOCK-LIFSHITZ SOLUTIONS 



In Ref . [15| , the authors define a finite action for Lifshitz theory in 4-dimensional Einstein 
gravity which satisfies 61 = with appropriate boundary conditions by adding appropriate 
local counterterms to the action (jlj). In this section we generalize this action to the case 
of (n + l)-dimensional Lovelock gravity. As in the case of Einstein gravity, to preserve the 
diffeomorphism invariance of the action the counterterms rendering the action finite and 
yielding a well-defined variational principle should be covariant in the boundary fields. We 
consider I = Jbuik + -^bdy, where /bulk is given in Eq. (jl]) and Jbdy is the sum of the boundary 
terms which are needed to have a well-defined variational principle and the counterterms 
which guarantees the finiteness of the action. Jbdy, for the case of zero curvature boundary 
which is our interest, may be written as 

1 /" .„ / — 7- r n — 1 / , a(n 



+3a3 - ^) + i/(.4„ A°) I + (20) 

where the boundary dAi is the hypersurface at some constant r, ha/s is the induced metric, K 
is the trace of the extrinsic curvature, Ka/s = V(q,?t,^) of the boundary (where the unit vector 
n'^ is orthogonal to the boundary and outward-directed), a{n) and b{n) are two dimensionless 
constants depending on n, and J and P are the traces of [l6| 

= ^{2KK^,K; + K.sK^'K^p - 2K^,K^'Ksp - K^K^p) (21) 
Pa/3 = ^{[K^ -6K^K''^K^s + 8KK^sK^^K'^ -6K^sK^'K,^K^^ + S{K^5K''^Y]K^^ 

+ {12K^ - 12K,^K''f')K^^K^^Ks^ + 24K^^K^^KseK'^Kp^} (22) 

In Eq. fl2Ul) . Idem is a collection of terms involving derivatives of the boundary fields, which 
could involve both the curvature tensor constructed from the boundary metric and covariant 
derivatives of Aa, Since the boundary is fiat and the fields are constants for [TTl this term 
will not contribute to the on-shell value of the action for the pure Lifshitz solution or its first 



variation around the Lifshitz background and therefore we ignore it throughout the paper. 
The matter part of the action fl20|) is the same as the matter part of the action in Einstein 



gravity. Thus, as in Ref. [15|, an arbitrary function is added to the action which 

is due to the fact that on the boundary A^A" = —q^ is constant for Lifshitz solutions. 

The variation of the action about a solution of the equations of motion is just the boundary 
term, 

1 



51 



IGvr 



dM 



flU'- 



(23) 



1 



+/'(/l„/l")(2A„M" + Ao^ApSh'^'') - -f{A^A'')h^p6h 



a/3 



where U^f, = U^^> + U]^> + H^J, with 



n 



(1) 

a/3 



K, 



al3 



n 



I 



-K 



n(2) 

■'-'■a/3 



a\n, 

2a2 I 3Jq,s — JhaB H — —h 



n 



(3) 
a/3 



'a 13 

3a3 ( 5P„/3 



'-a/3 



b{n 



al3 



As in the case of Einstein gravity 15|, one has for the Lovelock-Lifshitz spacetime ( llip 

n^'F^JA'' = zqr'/l'6A\ 



z-1 



and it is a matter of calculation to show that if one chooses 



a{n) = -(n - l)(n - 2)(n - 3), 
6(n) = (n — l)...(n — 5), 



then one also has 



That is. 



Equation 



n 



(2) 



0, n 



n 



(3) 



0, 



(2) _ 2{Z - 1)q2 2 



0, n 



(3) 



3(2; - l)d3 2 



/5 



r 6ij. 



r 6ii. 



shows that, as in the case of Einstein gravity 



(24) 



151], there are variations 



involving 6hij and 6A^ that we need to cancel. Using the same argument as in 



15| , the finite 



action may be written as 



/ d'^xV-hiK + 2a2J + SasP 



[n 



/5 



-^V^^^) (25) 



It is remarkable that fixing a single coefficient suffices to cancel both divergences associated 
with 6h'^^ and 6A^. Note that when z = 1, then q = and these definitions reduce to the 
familiar AdS rules for third order Lovelock gravity 16|. If one defines 



S. 



-h 



16n 



zq 
21 



levr 



[n>^F^^ + za{-A^A^)-'/^A, 



(26) 
(27) 



then the general variation of the action is 

61 = j d''x{So,p6h''^ + SJA"). (28) 

In the background (ITT]) , we have Sa/s = 0, = due to cancellations between the different 
terms, and this action satisfies 61 = for arbitrary variations around (fTTl) . Thus, we have 
a finite on-shell action which defines a well-defined variational principle for our background 
spacetime. 

In the holographic renormalization programme for gauge-gravity duality in relativistic 
field theories, one computes the finite stress tensor after constructing a well-defined finite 



action, since it carries important physical information a 



30ut the dual field theory. This 
221]. For asymptotically Lifshitz 



job has been done for asymptotically AdS spacetimes |21 
spacetimes, the dual field theory is nonrelativistic, so it will not have a covariant relativistic 
stress tensor. However one can define a stress tensor complex ISj, consisting of the energy 
density £, energy fiux Si, momentum density Vi and spatial stress tensor Vij, satisfying the 
conservation equations 



dt£ + d,£' = 0, dtV, + diV. 



0, 



where 



and 



£ = 2S' 



S'A, 



= 2S\ - S'At 



-2S\ + S'Ai V' 



-2S\ + S^Ai. 



(29) 
(30) 
(31) 



where in the Lovelock case 5*0,^ and are given by fl26f27l) respectively. Consequently the 
stress tensor complex (l30l) and (|3T|) satisfies the conservation Eqs. (|29|) . 



IV. THE CONSTANT C, 



In this section, we want to calculate the constant Co, which conserved along the radial 
coordinate r. Since there is no exact Lovelock-Lifshitz solution (except under special cir- 
cumstances), we calculate it at the horizon and at infinity. We will use this to relate the 
constant that appears in the expansion at r = oo to the coefficients at the horizon. 



A. Co at the horizon 



We begin by first reviewing the results found in 12|] for the expansion near the horizon. 



Requiring that /(r) and g{r) go to zero linearly, that is 

/(r) = /i {(r - ro) + ^(r - rof + Ur - r^f + U{r - r,Y + •••}, 
g{r) = giir- Tq) + g2{r - r^f + g-i{r - r^f + g^ir - Tq)^ + 

h{r) = fl'^ [ho + hi{r - ro) + /i2(r - ro)' + h:,{r - r^f + h^{r - r^Y + ...}, (32) 

and inserting these expansions into the equations of motion arising from the action (jlj) for 
the metric (IE]) with the conditions (JTOll [Eqs (2.8)-(2.10) of Ref. 12|], and solving for the 
various coefficients, we find that = 0. This is consistent with the fact that the fiux dA 
should go to a constant at the horizon. 

By rescaling the time coordinate, we can adjust the constant /i by an overall multiplica- 

1 /2 

tive factor (note the use of fi in the expansion of h{r) as well, which is due to dt in the 
one-form A). It is straightforward to find all the coefficients in terms of the two constants 
ro and hi. 



Using the above expansion, we evaluate the constant Cq ( 1181) at r = ro 



^ rp^W fi9i /oq\ 
L-o • (^^j 



This must be preserved along the fiow in r. 
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B. Co at infinity 



We now turn to the calculation of Co at large r. In order to do this, we investigate the 
behavior of the metric functions at large r by using straightforward perturbation theory: 

fir) = l + £/i(r), 
g{r) = l + egi{r), 
h{r) = 1 + ehi{r), 

and finding the field equations up to the first order in e. We obtain 

= 2r^h'l + 2{n + z)rh\ + zr {g[ - f[) + 2{n - l)zgi, 

= 2(z - l)rh\ + {n- l)rg[ + [z{z - 1) + n{n - 1)] g^ - {z - l){n + z - - 2h{), 
= 2(z - l)rh[ + (n - l)rg[ + [z{z - 1) + n{n - 1) + 2{n - l){z - l)B] gi 

-{z-l){z-n + l)Ui-2h) (34) 

where B = {l'^ — 4d2^^ + Qds) / L^. Note that all the parameters of Lovelock gravity are in B, 
with i3 = 1 in Einstein gravity. 
The solution of Eqs. is 

hi{r) 
Mr) 





C2 


C3 


;^n+2— 1 


^(n+2-l+7)/2 


^(n+2-l-7)/2' 




C2-P2 




^n+z— 1 


^(n+z-l+7)/2 


^(n+z-l-7)/2' 




C*2G2 


CaCs 




^(n+^-l+7)/2 


^{n+z-l-7)/2' 
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where 

7 = {(17 - 8B)z^ - 2{3n + 9 - 8B)z + + 6n + 1 - SB}^^^ , 
Fi = 2{z-l){z-n + l)}C'\ 
F2 = {Ti-T2){8zlC[{z-l)B + 2n + z-3]y\ 
F3 = {Ti + 72) {{z -l)B + 2n + z- 3]}"' , 
Gi = 2{z-l){n + z-l)}C-\ 
G2 = (Gi + Q2) {8zlC [{z -l)B + 2n + z- 3]}"' , 
G3 = {Gi - G2) [{z -l)B + 2n + z- 3]}'' , 
/C = {z-l){n + z-l)B + z{z-l) + n{n-l), 
J^i = 8{z - l)[{z - l){n + z~l)B + z{z - 1) + ^(n - 1)] 

X [{z - l)(n + 3^ - 3)B - + (n + 3)z + n{n - 2) - 1], 
J-2 = ^[n - 1 + {z - l)B]{^8il + B)iz - If 

+ {17n -9 + 8B){z- if + 2{n + 8)(n - 1)(^ - 1) + ^^(n - 1) - (n - 1)7 
^1 = 8(^ - 1)[2(2 - 1)B -3z + 3n- l][{z - l)(n + ^ - l)i3 + z(z - 1) + n{n - 1)] 
= |8(l + i3)(^-l)3 + (17n-9 + 8i3)(^-l)2 

+2(n + 8)(n - l){z - 1) + n2(n - 1)|7 - (n - 1)7^ 



The case of 2; = n — 1 in Einstein gravity needs special consideration. In this case 7 = 2z, 
and therefore the second term in Eqs. fl35|) is exactly the same as the first term. Also, one 
should choose C3 = in order to have suitable asymptotically behavior for the functions. 
It is a matter of calculation to show that the solution of Eqs. flMl) for z = n — 1, together 
with the constraint that these functions should go to zero for large r, may be written as 

"-^y' ) ^2(n-l) ' 

, , , 3n-A C2 



(r2-l)(2r2-3)r2("-i)' 
^ 2{n-2){G, + C2\nr) _ {n^ - 2)^2 

; - 3)r2(«-i) {n - l){2n - 3)2r2("-i) ' ^ ' 

where Ci and C2 are integration constants. 

Now we want to calculate the conserved quantity Cq. We restrict ourselves to the case 
that 7 > n + 2; — 1, which for Einstein gravity means that z > n — 1. However in Lovelock 
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gravity, this condition holds provided z> {n — B)/{2 — B), where B depends on the Lovelock 
coefficients. Note that this condition does not hold for B > 2. If z > {n — B)/{2 — B) then 
C3 = 0, since at large r the functions gi{r) and hi{r) should go to zero as r goes to 

infinity, and the contribution in Cq from C2 is zero. Thus, one has only the first terms in 
the expansions (153]) and the constant Cq can be obtained as 

^ 2{z-l){z + n- -n + 1)1' + 2{n - 2)a2l^ _ 3(n + z - 3)d3}^ 



We also find that the terms due to the second order perturbation do not contribute to the 
conserved quantity (ITTIl . 

The above constant (1371) in Einstein gravity for z = n — 1 is zero. In this case we must 
use the expansion (|36|) to obtain the conserved quantity Co, which gives 

4fn — 1) 

Co = (2n-3)/" ^' iz=n-l, Einstein) (38) 
V. BLACK BRANE THERMODYNAMICS 



The entropy of a black hole in Lovelock gravity is 23 1 

S = -^Y.^c,^ d'^-^^^gCk-u (39) 
fc=i 

where the integration is done on the (n — l)-dimensional spacelike hypersurface of the Killing 
horizon with induced metric g^^ (whose determinant is g), and Ck is the fcth order Lovelock 
Lagrangian of g^^. The entropy of a black brane per unit volume of the horizon in third 
order Lovelock gravity is 

S=\rV- (40) 
The temperature of the event horizon can be obtained by using the expansion (l32l) . 

yielding 

T = ^o^'v^ (41) 

Using the value of Cq at horizon ( 133|) and Eqs. ( HOj) and ( 14T|) . we find 

Co = IGnTS. (42) 

Alternatively, one can calculate the energy density of the black brane by using Eqs. (!26l) . 
([27D and ([30]). We find that 

(^n-l){z-l)iz + n-l){{z-n + l)l' + 2in-2)a2p -3{n + z-3)as}^ 
%^^zl^ ^^^^ 
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For z = n — 1 the above expression for the energy density in Einstein gravity vanishes and 
we must use the expansion f l36|) to find the energy density in Einstein gravity with z = n — 1. 
This gives 

Dividing Eq. (03]) by Eq. ([37]) (or, ior z = n - 1, dHD by ([38])) and using Eq. ^ we 
find 

T) — 1 

£ = — — —TS. (45) 

n+z-1 ^ ^ 

which is vahd for all z. Note that both the constant Cq and the energy density £ are different 
from their counterparts in Einstein gravity since the coefficient of r"*^""^^"^) at large r for 
the functions f{r), g{r) and h{r) differs. However, as one can see from Eq. (l45l) . £{T, S) is 
the same in both Einstein and Lovelock gravity. 

d equations have exact solution, the energy density. 



For the case of z = 1, for which the fie 
entropy density and temperature are 16| 



(n — l)m [n — l)rQ 

O — 



T = ^, s = Kr' 



AttP' ^ 4'° ' 

and therefore £ = [n — l)ST/n, consistent with Eq. (H5l) for z = 1. 



Using the first law of thermodynamics d£ = TdS with the relation (1^5!) for the energy 
density, one obtains 

logT = ^logS + r, (46) 

where F an the integration constant that depends on the Lovelock coefficients, z and the 
dimension of spacetime. It can be found numerically. This result is consistent with our 



previous numerical solution 12|], which shows that the slope of logT versus log 5* are the 
same for Einstein and Lovelock gravity, while F is different. Furthermore, by use of Eqs. 
lO]) and fll6]) . the temperature is proportional to Tq. 



VI. CONCLUSIONS 

In this paper, we introduced the finite action of third order Lovelock gravity in the pres- 
ence of a massive vector field with a fiat boundary. Indeed, we generalized the counterterm 



method introduced in (l6| for asymptotically AdS black branes to the case of asymptotic 
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Lifshitz black branes. We also defined the finite stress tensor complex, and computed the 
energy density of the Lovelock-Lifshitz black branes. We then used the field equations to 
find a conserved quantity alongthe r coordinate. This constant, which is the generalization 
of the constant introduced in 9[, has the role of connecting the metric parameters at the 
horizon and at infinity. We used these generalizations to investigate the thermodynamics of 
Lovelock-Lifshitz black branes introduced in We found that the relationship between 
the energy density, temperature, and entropy density is unchanged from Einsteinian gravity, 
even though the subleading large-r behavior of Lovelock-Lifshitz black branes is different 
from the Einsteinian Lifshitz solutions. We made use of the first law of thermodynamics to 
obtain the relationship between entropy and temperature. We found this is also the same 
as the Einsteinian case, apart from a constant of integration that depends on the param- 
eters in the Lovelock action. These results are consistent with the numerical analysis of 
Lovelock-Lifshitz black branes investigated in jl2|. 

Our counterterm method can be applied only to the case of A; = solutions. It would 
be interesting to generalize this method to the k = ±1 cases. Its generalization to qua- 
sitopological gravity with cubic-curvature terms that are not supersymmetric and therefore 
is different from the third order Lovelock gravity in the context of holography 14] would 
also be of interest. It is known that the third order Lovelock term does not contribute to 
the three-point functions of gravitons in fiat space j24|, while the corresponding contribution 



in the AdS background does not vanish 13|. Employing our action to calculate two and 



three-point functions in the pure Lifshitz background, and considering their differences with 
respect to the case of asymptotic flat and AdS solutions is another topic worth investigating. 
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